Analytical Model of Spin-Polarized Semiconductor Lasers 
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We formulate an analytical model for vertical-cavity surface-emitting lasers (VCSELs) with injec- 
tion (pump) of spin-polarized electrons. Our results for two different modes of carrier recombination 
allow for a systematic analysis of the operational regimes of the spin- VCSELs. We demonstrate 
that threshold reduction by electrically-pumped spin-polarized carriers can be larger than previ- 
ously assumed possible. Near the threshold, such VCSELs can act as effective non-linear filters of 
circularly-polarized light, owing to their spin-dependent gain. 



Spin-dependent properties have been successfully used 
in metallic magnetic multilayers for a variety of device 
applications exploiting magnetoresistive effects^ Unfor- 
tunately, much less is known about practical paths to 
implement spin-controlled devices that would go be- 
yond magnetoresistance. 1 ' 2 - One such encouraging de- 
velopment is demonstrated with spin-polarized vertical 
cavity surface-emitting lasers (spin- VCSELs).—! 4 - Spin- 
polarized carriers created by circularly-polarized photo- 
excitation or electrical injection can enhance the per- 
formance of spin- VCSELs as compared to their conven- 
tional (spin-unpolarized) counterparts. This work- 3 -* 4 - has 
demonstrated threshold current reduction and indepen- 
dent modulation of optical polarization and intensity. 
While numerical results using the rate equation (RE) 
description of spin- VCSELs were already presented,— i 4 - a 
large number of materials parameters makes it difficult 
to systematically elucidate how the spin-dependent ef- 
fects modify the device operation. Consequently, even 
the widely accepted theoretical limit of maximum 50% 
threshold current reduction needs to be re-examined. 

For clarity of our approach we formulate a simple RE 5 
model of spin- VCSELs which, while similar to prior nu- 
merical work; 3 -! can be solved analytically. We consider 
a quantum well (QW) as the active region in the VC- 
SEL. Spin- resolved electron and hole densities are n±, 
P±j where +(— ) denotes the spin up (down) component; 
the total carrier densities are n = n + + n_ , p — p + + p_ . 
Analogously, for photon density we write S — S + + S~ , 
where +(— ) is the right (left) circularly-polarized com- 
ponent.^ Electrically or optically injected/pumped spin- 
polarized electrons into the QW can be represented by 
current density J = J + + J_ (normalized to unit charge 
and an effective volume in the QW region 4 -). Typically, 
the spin relaxation time of holes is much shorter than 
for electrons ?i rf <C r™, implying that the holes can be 
considered unpolarized with p± = p/2. 

In conventional VCSELs, the optical gain term, de- 
scribing stimulated emission, can be simply modeled as 
g(n, S) = go(n — n tlan )/(l + eS), where go is the density- 
independent coefficient f 4 - ntran is the transparency den- 
sity, and e is the gain compression factor. We general- 
ize this relation in the spin-polarized case as g(n, S) — * 
9±{n±,p±,S ± ) = go(n±+p±-n tian )/(l+eS). This form 
differs from the previously employed gain expressions, 3,4 



as it explicitly contains the hole density, but it coincides 
with the rigorously derived gain expression from semi- 
conductor Bloch equations^ 

The charge neutrality p± = p/2 = n/2, allows us to 
recover the spin-unpolarized limit for the gain expression, 
as well as to decouple the REs for electrons from those for 
holes. The spin-polarized REs for electrons thus become 

dn±/dt = J± - g±(n±, S^)S^ - (n± - n T ) /r s n - J^,(l) 
dS ± /dt = Tg T {n Tl S ± )S ± - S ± /r ph + f3TRf pl (2) 

where r p h is the photon lifetime, T is the optical con- 
finement coefficient, (3 is the spontaneous-emission cou- 
pling coefficient, and R^ p is the radiative spontaneous 
recombination rate. For the spin-unpolarized case, R^p 
is typically assumed to be quadratic in carrier densi- 
ties,— R sp = Bnp/2, where B is a temperature-dependent 
constant^ Another simple form, R sp — n/r r , where T r 
is the recombination time, is a good approximation at 
high n,— but was not considered in the prior work on 
spin- VCSELs. 3,4 For spin-polarized electrons these forms 
are generalized as^S Rf p — Bn±p± = Bn±n/2 and 
R^p = n±/r r . We refer to them as quadratic and linear 
recombination (QR and LR), which can be viewed as the 
opposite limits of a general expression for R^ p (n±,p±).— 



To describe the solutions of Eqs. (|1|2|) . which are the 
basis of our theoretical approach, we introduce polar- 
izations of injected electron current Pj = (J+ — J_)/J 
as well as of electron and photon densities P n — (n + — 
n_)/n, Ps = (S + — S~)/S. Electrical injection in QWs, 
using Fe or FeCo, allows for^ 2 \P n \ - 0.3 - 0.7 with 
similar values for \Pj\, while \P n \ — > 1 is attainable op- 
tically at room temperature.- 1 ^ 3 - In the unpolarized limit 
J+ = J_, n+ = n-, S+ = S~ or P., = P n = P s = 0. 
As in the recent experiments? 4 - we consider a continuous 
wave operation of VCSEL and look for steady-state solu- 
tions of Eqs. (|1)2|L Guided by the experimental ranged 
of (3 ~ 10~ 5 — 10 -3 , we mostly focus on the limit = 0, 
for which all the operating regimes of spin- VCSELs can 
be simply described, and, additionally, consider e = 0, 
relevant for moderate pumping intensities. 

We first show analytical results for unpolarized [Pj = 
0) VCSEL and LR in the inset of Fig. [1] Injection current 
density is normalized to the unpolarized threshold value, 
Jt = Nt/t t , with Nt denoting the total electron den- 
sity at (and above) the threshold, Nt — n(J > Jt) = 
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FIG. 1: Photon densities of the left- (S~) and right- cir- 
cularly polarized light as a function of electron current J with 
polarization Pj = 0.5, infinite electron spin relaxation time 
and linear recombination. J is normalized to the unpolarized 
threshold current Jt and to St = Jt IVph. The vertical 
lines indicate Jti and Jt2 thresholds. Inset: total photon 
density (S = S + + S~) is shown for spin-unpolarized laser 
(Pj = 0) and two spontaneous-emission coupling coefficients 
((3 = 0, 0.002), as well as a spin laser with Pj = 0.5, (3 = 0. 



(TgTph) 1 + fitran, while photon density is normalized to 
St = Jr^Tph- A small difference in S(J) between the 
vanishing (solid line) and finite (3 (dotted line, f3 = 0.002 
overestimates the experimental values^) shows the ac- 
curacy of (3 = approximation. The unpolarized case 
has two regimes: for J < Jt the device behaves as a 
light-emitting diode (LED), with negligible stimulated 
emission; for J > Jt, it is a fully lasing VCSEL. 

With finite Pj and r™ — > oo, we reveal a more com- 
plicated behavior further explored in Fig. [TJ main panel. 
The two threshold currents, Jti and Jti {Jti < Jt < 
Jt2, the equalities hold only when Pj — 0), delimit three 
regimes of a spin- VCSEL. (i) For J < Jti it operates as 
a spin-LED^ii^ (ii) For Jti < Jt < Jt2, there is mixed 
operation: lasing only with left-circularly polarized light 
S~ (we assume^ J+ > J-), which can be deduced from 
the spin-dependent gain term in Eqs. (jl!2p . while S + is 
still in a spin-LED regime (S + — * for (3 — > 0). (iii) For 
J > it is fully lasing with both S ± > 0. 

RE description of spin-unpolarized lasers reveals that 
the carrier densities are clamped above Jt- We have 
found a related effect in the spin-polarized case. For J > 
Jti , Eq. @ for S~ can be divided by S~ , showing that 
the quantity n + + p+ = n+ + n/2 (we recall the charge 
neutrality condition and t£ — * 0) will be clamped at Nt- 
If Jti < J < Jt2 then neither n + nor rt_ are separately 
clamped. If J > Jt2, then n_ + n/2 = Nt must hold 
in addition to the previous condition n + + n/2 = Nt- 
These two conditions yield n± — Nt/2, independent of 
Pj. Thus, above Jt2, the sum of Eqs. @ for S~ and 
S + reduces to the usual unpolarized equation, and S is 



independent of Pj (inset of Fig. [TJ. However, if Pj ^ 
then we still find S~ ^ S + (Fig. [TJ J > J T ). 

Most of the results discussed above do not change qual- 
itatively for a finite r™ or QR. The general features, such 
as the existence of three regimes of operation, remain the 
same. However, the photon and carrier densities (for a 
given J) as well as the threshold Jti depend quantita- 
tively on the spin-flip rate and the recombination form. 
We investigate this dependence in Fig.r2]which shows the 
evolution of photon and carrier polarizations with the in- 
jection current. We consider both LR and QR forms and 
express our results using the ratio of the radiative recom- 
bination and spin relaxation times: t = t t /t™ . For QR 
the unpolarized threshold is Jt = BNjl/2, and we define 
r r = Nt/Jt by analogy with the LR case. 

When < J < Jti, then P n (J) depends on the re- 
combination form; P n is constant for LR, while for QR 
it grows monotonicallyJ^ Only for a very long spin re- 
laxation time, t — > 0, P n (J < Jti) — > Pj- At = 
there is no stimulated emission, so that Ps = 0. When 
Jti < J < Jt2, Ps = — 1 f° r an Y t, recall Fig. [TJ 
With only a partially polarized electron current, the spin- 
VCSEL emits fully circularly-polarized light, analogous 
to the spin- filtering effect in magnetic materials. In the 
same regime, P n decreases with J because, while n_ 
grows, n+ must drop in order to maintain n+ + n/2 = 
3/2n+ + l/2n_ = Nt- 

For Pj ^ 0, the Jt2 — Jti interval is widest for 
t™ — > oo. When r™ decreases then the lower threshold 
Jti grows towards its unpolarized counterpart Jt, and 
Jt2 — Jti contracts, Fig. [2j We find J T 2 = Jt/(1 - Pj), 
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FIG. 2: Electron (P n ) and photon (Ps) polarizations as a 
function of injection current J, with polarization Pj = 0.5 and 
different ratios of recombination and electron spin relaxation 
time, t. Linear (LR) and quadratic (QR) recombination are 
shown. J is normalized to the unpolarized threshold Jt- 

independent of r™, this should hold up to very small 
values of r™ (t ~ 10 2 ). In the extreme case of even 
smaller r™, for Jti < J < Jt2 the majority and mi- 
nority spin densities become almost equal; n + > Nt and 
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n_ < Nt, so that f3 ~ 10~ 4 can effectively drive the 
device to the fully lasing regime. For typical values 3,4 
of t ~ 5, setting j3 — is an accurate assumption. At 
Jt2, we recall n + = n_(= Nt/2) so that the spin- flip 
term in Eq. ([1]) for n_ is zero. The equation reduces to 
J- = 1/2(1 — Pj)J = R~ p , explaining why Jti depends 
only on Pj, but not on t. For the same reason, Pg and P„ 
are independent of t when J > Jti ■ The quantity Ps in- 
creases (\Ps\ decreases) with J as Ps — —PjJ/ ( J— Jt) to 
the asymptotic value Ps = — Pj. Thus Pj can be inferred 
from Pg only at sufficiently high injection. In all three 
regimes, defined by Fig. [TJ we find that < P n < Pj, 
(only for r™ — > 00 and J < Jti, P n — > Pj). For J > J7-2 
we find P n = 0. Therefore P n is not enhanced in the 
spin-VCSEL's active region. 

Experiments have demonstrated 3 -! 4 - that injecting spin- 
polarized carriers reduces the threshold current in a VC- 
SEL, i.e., Jti{Pj ^ 0) < Jt- We quantify this threshold- 
current reduction 3 -! 4 - with R = (Jt — Jti) I ' Jt, obtained 
analytically from REs 

LR: R = |P/ 1 /(2 + \P,j \ + 4i), (3) 
QR:P=l-2/(2+|P 7 |) 2 x (4) 

1 + 2 |P/| t + 4t 2 + (1 - 2f) y / (l + 2t) 2 + 4|P J |i , 

and we find R > 0, for any Pj 5^ 0, r™ > and both 
LR and QR. As expected, R — > for r™ — > (i — > 00), 
independent of the recombination form, since we have 
always assumed = 0. We illustrate further R(Pj) for 
both LR and QR in Fig. [3] We choose t — 1,5, close to 
the typical values r r ~ 1 ns and r s ~ 100 ps for QWs 
at room temperature^ From Fig. [3] or Eqs. (|3I4[) one can 
see that R(Pj,t) for LR is always smaller than for QR. 
For \Pj\ = 1 and r ™ ->■ 00 (but still rf = 0) we obtain 
the maximum current reductions 

LR: R = 1/3, QR: R = 5/9. (5) 

Remarkably, for QR the maximum R is larger than previ- 
ously assumed possible R = 1/2^ even though holes are 
completely unpolarized (rf = 0). This can be explained 
as follows. First we calculate n from Eq. Q for S~. 
The threshold density is reached when n + + n/2 = Nt, 
which gives n = (2/3)Nt for Pj — 1, i.e., for n + = n. 
Thus the threshold electron density is reduced by only 
1/3. Assuming S~ = at the threshold, the threshold 
current from Eq. ((T|) is given by J + = Bn + n/2. If Pj = 1 
then J = J + — Bn 2 /2 = (2/3) 2 J T , which yields the 5/9 
current reduction, a direct consequence of the quadratic 
dependence of recombination on n (note that for LR both 
the Nt and Jt are reduced only by 1/3). 

For comparison, we calculate R for QR in a priori 
"ideal" case, where n± — p± for any injection J and 



both t™,t% — > 00. Therefore the gain and QR terms in 
Eq. (U) for n + (n_) and in Eq. © for S~ (5+) must 
be modified by replacing n/2 with n+ (n_). Surpris- 
ingly, we obtain P(|Pj| = 1) = 1/2, a smaller reduction 




injection polarization \Pj \ 

FIG. 3: Threshold current reduction R as a function of in- 
jection polarization Pj. Shown are both linear and quadratic 
recombination (LR and QR) for various t, the ratio of recom- 
bination time and electron spin relaxation time. The hole 
spin relaxation time rf = 0, except for the dashed line and 
what was previously assumed "ideal" case, in which rf — > 00. 



than in the case of unpolarized holes. A simple calcula- 
tion yields n = Nt/2, i.e., the reduction of the thresh- 
old electron density is larger than for rf = 0. However, 
the interband recombination P+ is more efficient, be- 
cause none of the holes undergo spin-flip. Thus we obtain 
J + = J T1 = Bn\ = BN T /A = J T /2. For rf = and 
Pj = 1, half of the holes are inaccessible for recombina- 
tion with the fully spin-polarized electrons. This lowers 
the injection necessary to overcome the recombination 
losses and the interplay of stimulated and spontaneous 
recombination leads to a smaller R in the "ideal" case. 

In this work we show that even a simple rate equation 
model can reveal surprising trends for operation of spin- 
VCSELs. The maximum threshold current reduction is 
not achieved for infinite electron and hole spin relaxation 
times, but rather when the hole spin lifetime vanishes. 
The corresponding threshold reduction exceeds the pre- 
viously considered theoretical limit. We expect that the 
transparency of our analytical approach will help to elu- 
cidate the operation of spin-controlled lasers, which con- 
tinue to be actively studied. 15 
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